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1 Introduction 


Looking at the great success that the Standard Model has had since it was 
introduced more than 20 years ago, it is quite striking that we still lack under¬ 
standing of the strong interaction part in the low energy regime. In spite of 
excellent results in the perturbative QCD we are unable to produce any ana¬ 
lytical computation of quantities such as the magnetic moment of proton which 
is known with great accuracy. Some important puzzles, like where does Aqcd 
come from, still need to be addressed. Apart from the numerous lattice results 
(which reinforce our belief that QCD is the right theory for the strong interac¬ 
tions) the situation remains largely the same more than 20 years later. Many 
attempts to apply a variational approach have not yet produced any effective 
calculational method for the solution of the problems mentioned before. 

It was suggested in [Q that perhaps we need an alternative way to look at the 
Yang-Mills theory based on a more geometrical point of view. Namely, one tries 
to study quantum mechanics of the fields in the space of gauge-inequivalent con¬ 
figurations. The following analogy with quantum mechanics is used. Consider 
a free particle inside a box of size L. The lowest eigenvalue of the Hamilto¬ 
nian (which is just a Laplacian) is of the order of ^ This is realized by the 
state of the longest possible wavelength, X ^ L. It is clear therefore, that the 
spectrum of such a system will have a gap due to the fact that L is finite. In 
other words, the spectrum is going to be discrete as long as distances in the 
configuration space cannot become arbitrarily large. Feynman’s suggestion was 
to try to adapt this idea to Yang-Mills theory. 

In general, in order to determine the distance between field configurations 
we need to know the metric of the configuration space. The geometry of the 
configuration space of Non-Abelian gauge theories was considered by many au¬ 
thors (|^ 1^, Q) and, recently, in |^. The major problem is to extract the metric 
on a space of the gauge-inequivalent configurations A/G- 

A natural distance between two arbitrary gauge configurations Ai and A 2 
(that is gauge connections modulo gauge transformations) is given by 

\\A,-A2\f=znf{[TriiA{{x)-A2{x)f% ( 1 ) 

J X 

where Af is a gauge-transformed Ai. However the extremum solution of this 
expression is highly non-local and very difficult to work with. Nevertheless one 
can try to see some qualitative features of this distance. In ||^ it was argued 
that due to the non-Abelian nature of the gauge group the low-potential-energy 
configurations are in a space of finite diameter for 2-1-1 dimensional Yang-Mills 
theory. According to our very crude analogy with the quantum mechanics it 
would mean that the kinetic energy operator will have a discrete spectrum. Of 
course, these are only qualitative arguments that require a rigorous mathemat¬ 
ical formulation if one is to draw any conclusion on the existence and value of 
the mass gap. Also, as a word of caution we must say that some properties of 
the configuration space obtained in 2-1-1 dimensions may not necessarily be the 
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same as in the 3+1 dimensional case (as, for example, the statement above in 
the case of YM in 3+1 and 1+1 dimensional sigma model |^, ^). 

Another reason to believe that this approach could lead us to a better un¬ 
derstanding is the recent progress in Hamiltonian formulation of the 2+1 di¬ 
mensional Yang-Mills theory. In a series of papers 01 it was shown that by 
introducing special gauge-invariant variables one can prove that the (properly 
regularized) volume of the configuration space for the non-Abelian theory is 
finite, while the corresponding quantity for the Abelian field is infinite. The 
discreteness of the spectrum of the kinetic term, has been also shown ex¬ 
plicitly and the string tension computed in is in remarkable agreement with 
the recent Monte-Carlo simulations |^.) 

Behind all this discussion a natural question arises: What is the behavior of 
the configuration space as we integrate out high momenta degrees of freedom? 
Of course we can not answer completely this question, and even a partial re¬ 
sponse deserves a profound, probably non-perturbative, analysis. However we 
still can say something, and answer specihc matters as the behavior of the con¬ 
figuration space metric under the renormalization group, which encodes much 
of the properties of the space. Notice that those issues acquire a special signifi¬ 
cance in Hamiltonian formalism: there the kinetic term is essentially a Laplacian 
in configuration space and the energy spectrum, at least in the strong coupling 
regime, is dominated by it. 

The choice of this topic was not fortuitous but we have an idea in mind. In 
the standard understanding of asymptotic freedom all the signihcance is put in 
the interaction potential through the statement that the couplings decrease to 
zero as the energy at which the theory is tested increases to infinity. However 
in the Hamiltonian picture we can state the problem in a different way. In the 
Hamiltonian there is an obvious competition between the kinetic energy and the 
potential energy. So it is natural to analyze the asymptotic freedom in its “dual” 
form, i.e., the variation of the kinetic energy (with respect to potential energy) 
as the renormalization scale changes. In particular we conjectured the follow¬ 
ing, in asymptotically free theories the effective distance between configurations 
decreases as the momentum is lowered. 

Notice that this question does only make sense in a Hamiltonian formalism, 
since it is only in this case that all the geometrical properties of the config¬ 
uration space can be precisely defined (the Hamiltonian is an operator that 
acts on the Hilbert space of functionals defined on the configuration space). 
It can be argued that one can always work in a Lagrangian formalism and at 
the end construct the respective Hamiltonian by the standard Legendre trans¬ 
formation. However this cannot be done if renormalization is involved at any 
moment. In the Lagrangian formalism, even in the framework of perturbation 
theory, the process of renormalization generates high time-derivative terms in 
the Lagrangian, making impossible even the very definition of a Hamiltonian, 
as the system is governed by high order time-derivatives equations of motion. 
Even an application of the Ostrogradsky method can in general lead to negative 
norm states if higher time derivatives are involved in the Lagrangian. Thus, a 
prescription for renormalization within the Hamiltonian formalism is indispens- 
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able. One must not leave the phase space, thereby maintaining the first order 
time derivative nature of the equations of the system. 

Therefore we will use a novel procedure of renormalization of Hamiltoni¬ 
ans, already introduced in ref. | p^ where a successful renormalization of the 
Hamiltonian of A was performed. This method relies on the successive diag- 
onalization of the Hamiltonian by performing iterative unitary transformations 
and subsequent projection onto the Hilbert space of low energy modes. It re¬ 
sembles, in spirit, the renormalization approaches of Glazek and Wilson |pd| , 
Wegner El and several other authors j^, |^, [l^ |^, though in practice ap¬ 
pears very different. To support our conjecture we have done a detailed analysis 
of two particular examples: Quantum Electrodynamics and Yang-Mills theory 
in 3-1-1 dimensions, where we constructed the renormalized Hamiltonian up to 
one loop. 

The paper is organized as follows: In Section 2 we present briefly our renor¬ 
malization group technique for Hamiltonians. In Sections 3, 4 and 5 we compute 
the 1-loop renormalization of the SU{N) Yang-Mills Hamiltonian. In Section 6 
we did the same with Quantum Electrodynamics. In Section 7 we discuss the 
relation between renormalization of the metric and asymptotic freedom. Finally, 
Section 8 contains some conclusions. 


2 Renormalization in the Hamiltonian formal¬ 
ism 

Let us review in this section the formalism of Hamiltonian renormalization in¬ 
troduced in ref. [ p^ . 

Consider a system described by Hamiltonian H which has already some 
large cut-off A built into it; that is the system is ultraviolet-finite from the 
very beginning. We are assuming that the Hamiltonian is written in terms of 
renormalized fields and couplings up to the scale A and incorporates all the 
renormalization Z{A) factors. 

Now if fj, is some intermediate scale, /r < A, we want to find the effective 
Hamiltonian H(fj.) that has the same low energy spectrum as the original Hamil¬ 
tonian. More precisely, we want to find the operator that has the same spectrum 
as the original Hamiltonian when projected onto the Hilbert subspace generated 
by the excitations with frequency less than /r: 

Hitj.)=Vio^H{A)Vio^. ( 2 ) 

where Viow is the projector onto the low frequency subspace. 

We will show that, in the framework of perturbation theory, it is possible to 
partially diagonalize the Hamiltonian H (A) and construct the vacuum state for 
the high frequency modes. Thus, the low energy effective Hamiltonian would 
take the form 


= {0h^gh\UHA,^i) H{A) U{A,fi)\0h^gh)■ (3) 
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Notice that it is not necessary to diagonalize completely the Hamiltonian as 
we only need to identify the low energy states and not the whole spectrum. 

Suppose that after a unitary transformation we bring a Hamiltonian to the 
form 

H(low) + Hfree(high) + V (4) 

where H(low) contains only low frequency operators, Hfree(high) is a free Hamil¬ 
tonian for the high frequency modes and V has the special form 

y = ^ al Skp{fi,a\a) ap (5) 

k,p>^ 


with SkpifJ., a^, a) an arbitrary operator of low and high frequency modes. Thus, 
using standard Rayleigh-Schroedinger perturbation theory it follows that the 
correction to an arbitrary state |n) is given by 


|fc) = 

l^n 


(^1^1^) in , {l\V\m){m\V\n) 


( 6 ) 


Therefore, if |n) is the vacuum state of TJfree(high), it is annihilated by V and 
there is no correction to it at any order in perturbation theory. 

Now we will show how to find the unitary transformation U{p, A) that trans¬ 
forms the Hamiltonian into the form (^. First we split the Hamiltonian in four 
pieces: 

H = Hi+H2 + Va + Vb. ( 7 ) 

Here Eli contains only the modes with k < fj,^ iJ 2 is the free part for all the 
modes with k > fj.^ Va contains the “pure” terms that have only high frequency 
creation operators or high frequency annihilation operators, but not both, and 
Vb the remaining terms containing at least one annihilation and one creation 
operator of high energy modes (we assume here that Va and Vb are normal- 
ordered with respect to the free high frequency vacuum). Notice that the term 
Vb is already in the form (||). Then, according to the preceding discussion, the 
objective of the renormalization procedure is to “eliminate” the terms in Va 
that spoils the form d,!)- 

We will proceed iteratively: we write the unitary operator in the form U = 
U 0 U 1 U 2 ■ ■ •, and compute each Un perturbatively. Each factor Un will bring the 
Hamiltonian to the form (^) at a given order in /i/A. Therefore, there will be 
two expansion parameters in this procedure: one is the coupling constant of the 
theory A and the other is the ratio /i/A (effectively oJiow/'^high)- 

Let us parametrize t/„ as exp{i^ln), where is a hermitian operator to be 
determined. We start from equation (^) and perform a first unitary transfor¬ 
mation expanding in powers of H (in the general case H is at least of order A so 
at a given order only a finite number of terms are needed): 

g-zOo ^h2VVaE Vb)E^° = Hi+H2 + VA + VB+i [i7i, Hq] + 

-I- i [E[2, Hq] -I- i [Va; ^ 0 ] ■ ■ ■ (8) 
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We want to eliminate the term Va which is of the first order in A, so fio has 
to be of the same order. Furthermore we want to generate an expansion in /i/A, 
so we impose the following condition on Qq 

tlH 2 ,fiol + VA = 0. (9) 

This equation can be solved perturbatively and since commutators with H 2 
generate time derivatives we have the desired expansion. 

Notice that equation (||) defines flo up to the terms that commute with H 2 . 
As it is shown in the appendix A this freedom corresponds to the arbitrariness 
of the definition of the low energy Hamiltonian up to a unitary transformation. 
Therefore we will assume some kind of “minimal” scheme - namely that Hq, 
after normal ordering, does not have a part that commutes with H 2 (terms that 
are functions of a\ak)- 

After Hq is chosen to cancel V 12 in the effective Hamiltonian a new mix¬ 
ing term of order A has appeared from Notice however that this 

new term is of higher order in fj./A and will be eliminated by the next unitary 
transformation Ui = Explicitly, 

+H2 + VA + = H1+H2 + VB + 

i [i?i, Hq] -I- i [H 2 , Hi] -|- i [iJi, Hi] -I- second order terms • • • (10) 

and we now choose Hi so that 

^[i^l,Ho]+^[i^2,Hl] = 0. (11) 

Using equations (0) and (^^ we obtain: 

+H 2 + Ui2)e*^«e*^i = Hi+H2 + Vb + ^ [U4,Ho] + 

i [Vb, Hq] -I- i [Vb, Hi] -|- i [Hi, Hi] — - [[Hi, Hq] , Hq] — 

i [[Hi, Ho], Hi] — i [[Hi, Hi], Hi] -I- higher order terms • • • (12) 

The obvious next step is to introduce H 2 in order to cancel i [Hi, Hi] and con¬ 
tinue with the same process. Then a simple question emerges: Where should 
we stop? To answer this question we have to consider the divergence properties 
of the terms introduced by each new H into the Hej f. These contributions like, 
for example, — \ [[Hi, Hq] , Hq] may diverge as A ^ 00. Nevertheless, the degree 
of divergence of each new term will be smaller as we introduce more and more 
H„ factors. In general, the following power counting can be used: since H„ is 
determined from 

[H2,H„+i] = -[Hi,H„] (13) 

and H 2 ~ A, Hi k, then 

Hn-l-l ~ 
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Therefore the next fl will be less divergent and eventually all the new terms 
introduced by this prescription will be convergent starting from certain n. At 
this point we will stop since for our purposes we are only interested in divergent 
contributions. 

Of course so far we have only eliminated the high momentum degrees of 
freedom up to the first order in the coupling constant. Requiring the absence of 
the - order mixing terms will lead to the introduction of a whole new series 
of unitary transformations, and the same arguments can be applied to them. 

Finally, the actual process of renormalization is performed by choosing the 
renormalization Z{K) factors of the original Hamiltonian to cancel the diver¬ 
gent contributions coming from evaluation of (i- We now turn to the explicit 
computation for QED and QCD Hamiltonians. 


3 SU{N) Yang Mills theory: Preliminaries 


The kinetic energy term for the Hamiltonian of Yang-Mills theory is essentially a 
Laplace operator on the configuration space. However, due to the gauge invari¬ 
ance, the physical degrees of freedom belong to the space of gauge connections 
modulo the group of gauge transformations, i.e, the space of non-equivalent 
gauge potentials, and a satisfactory parametrization of this space is needed. In 
this section we find a perturbatively adequate coordinate system of the configu¬ 
ration space (a similar analysis was done in ref. |^) and compute the associated 
metric . 

We consider a SU{N) Yang-Mills theory in the temporal gauge, Aq = 0. The 
canonical variables are the vector potential Af{x) and the electric held E°‘'^(x). 
They satisfy canonical commutation relations 

[A^{x),E'^^{y)]=z5<^Hl5^^Kx-y) (15) 


which permit the representation of the electric held as 


E‘^\x) 


S 

*5A“(a;)' 


(16) 


The Hilbert space of the theory is supplemented by the Gauss law that enforces 
a constraint on the wave functionals, essentially only allowing gauge invariant 
conhgurations, 

G'^ixMA] = -^D[A]f = 0 (17) 

where — e/“^°A°(a;) is the covariant derivative. 

The gauge potential can be written also as a Lie algebra valued held Ai (x) = 
t°'Af(x) where are the hermitian generators of the Lie algebra of SU{N) in 
the fundamental representation, normalized to tr(t“t^) = ^5°"^. 

The Hamiltonian can be written as 



d^x {E°-\x)E‘^\x) + B°-\x)B°-\x)) 


(18) 
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where we used the magnetic field i?f: 

B-\x) = (^d,AUx) + ^er^^A’;ix)Alix)^ . ( 19 ) 

Note that due to the constraint ( 0 )> not all the degrees of freedom in the 
Hamiltonian (|^) are “physical”. In order to isolate the physical degrees of 
freedom we will perform a change of coordinates in such a way that the Gauss 
law takes its simplest form. We parametrize an arbitrary configuration in terms 
of a “gauge fixed” configuration plus gauge transformations. The former will 
define a coordinate system of the orbit space. We write 

Ai{x) = g~^{x)Ai{x)g{x) + ig~^{x)dig{x) (20) 

where 5 is a S'C/(iV)-valued matrix and Ai is a configuration satisfying the 
Coulomb gauge condition 

y ^ A ^ = 0. (21) 

It is well known that the parametrization (^) with condition (^) is not well 
defined globally due to the Gribov ambiguity problem. However, we will work in 
the framework of perturbation theory where the parametrization d^-l^ defines 
an acceptable isomorphism in the configuration space. 

From eq.(pO|) we deduce, 

6Ai = g~^ [5Ai + [D[A]i,i5gg~^]) g ( 22 ) 

where D[Ai] = di — ieAi- 

Thus we can decompose the metric of vector valued configurations in “gauge 
fixed” and “pure gauge” parts: 

6s^ = J d^x SA^{x) SA^{x) 

= J d^x {SA^x) SA^x) + D[A\i\i5gg-^fD[A]T{i5gg-y) .(23) 

The normalization of the wave functionals is then given by 

= vkc J 

= J[DAt] Vg «'2[^] (24) 

where G is the matrix metric defined by equation (^). 

Using the Coulomb gauge condition eq.(|^) we can invert eq.(p^ and write 

* {Sgg-^ = (v • V • {SAR-y (25) 

and 

6A^ = - D^^[A] (v • D[A]^ {6AjR-^)K (26) 
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where R = R{g) is the adjoint representation representative of g. 

Notice that not all the components of Ai are independent as they are sub¬ 
ject to condition so we can parametrize the space of gauge configurations 
modulo gauge transformations with the 2(N^ — 1) functions .4?, i = 1,2. 
From eqs.(p5[-p^ we can write the functional derivative as 


J R-'-- (s„S‘“S(x.v) + (D[A]‘Vf(3 . V)-;^,) 


K9~^Sg)^{y) 



+ 


(54' 

SA^) 


(27) 


where Pij(x, y) = 6ij6(x, y) — ViVj A ^(x, y) is the projector on the transverse 


modes. 

As previously announced, in these variables the Gauss law has the simple 
form 


J4'[A] 
{iSg g-^ f 


(28) 


so it is trivially imposed by demanding that the wave functionals be independent 
of g. 

It is useful to dehne a “transverse” functional derivative as 


,5^ 

SA^ix) 





z) 


6 

(5A|(z) 


(29) 


restoring rotational invariance at the expense of modifying the canonical com¬ 
mutation relations: 


^nx),Tinr 


SAAy) 


= -Pi]{x,y) S' 


ab 


(30) 


Now we can write the kinetic energy term in terms of the variables Af (taking 
into account eq.(p^)), 


(4'i|T|^2) = 


1 


Vol G 

J[DA'^ 


[DA'l] ; 

VdetG 


SAf{x) dAf{x) 
(fx Sy (fz X 

-1 


+ (d[X1JV-(5 . 

+ (D\A\V^t(D ■ V),-y,)“) 


1 ' 
(5^ 4-2 


SAUzY 


That is, the kinetic energy density term takes the form 

1 ^{ia)(jb) , 


T = - 


2VG 


£^‘^ix)VG G^X/^=^\x,y) S^\y) 


(31) 


(32) 
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cT 

where 5“ = and G'_ 4 /g is the effective metric of the space of gauge 

configurations module gauge transformations, 




+ [diay^vud ■ 


be 


(33) 


As we mentioned above, the previous analysis was only valid in the frame¬ 
work of perturbation theory where condition & defines a local system of coor¬ 
dinates on the orbit space. So it is consistent with this approach to compute all 
the elements of T (the metric of the space of gauge transformations G and the 
metric of the true configuration space G^/g) in powers of the coupling constant 
e. In fact, after a straightforward computation, the kinetic energy, up to order 
e^, can be shown to be 

T = ^Jd^xd^yG^X/^^^\x,y)£<^\x)£G{y)- 

A“^(0) J d^x Af£°‘'^{x) (34) 


where ca is the Casimir of G in the adjoint representation, 

G%^^^^\x,y) = 6^^S<^^6ix,y) - e^r‘^f^^A‘!ix)A^^{y)A-\x,y) + 0{e^) (35) 


and the transverse variables Af and £“ satisfy the commutation relations (see 


eq.(m): 


[Anx),£^{y)]=tS‘^^ P.,{x,y). 


(36) 


4 Renormalization of Yang-Mills Hamiltonian 

In this section we will compute the renormalization contribution to the and 
terms in the effective Hamiltonian. In what follows, we will assume the 
“gauge-fixed” variables from the previous section and will use A and E instead 
of A and £. Using expression (^^ the gauge-fixed Yang-Mills Hamiltonian can 
be written as follows: 



a^{x,y)Ef{x)E]{y) + \icAA‘^i{x)E'l{x)G{x,x 


+ -Bt{x)Bt{x) 
(37) 


where ca is N for SU{N) and 


o^ij{x,y) 

G{x,y) 


-y) + e^r'^^f‘^^A‘i{x)A]{y)G{x, y) + ■ ■ ■ 



-iq(x-y) ^ d 

( 2^)3 


(38) 

(39) 
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This is the so called “bare” Hamiltonian. For the loop calculations we have to 
introduce the Z-factors by the following procedure. Using let us rewrite 
expression as 

H = + ^ZAAEE{A)e^r'^^f‘^^AUx)A^^{y)Gix,y)Enx)E^iy) 

1 

+ -Zb 4A)B^ + eZi{A)r^^d,A"A\A‘^^ + — + • • • (40) 

Each of the Z-factors will have the following form 

Z=l + e/i(A) + eV 2 (A) + --- (41) 


The functions fn will be chosen order-by-order from the requirement that after 
integration of the modes from fj, to A all the corrections sum up in such a way 
that /n(A) ^ fnifJ-) and Z(A) ^ Z{fj,), accordingly. When doing the one- 
loop corrections one can therefore assume that all the Z’s are initially 1 and 
choose the corresponding f’s from the condition that high-cutoff dependence be 
cancelled after computing the i?e// to one loop. Since one-loop wave function 
renormalization in QCD is of the second order in coupling constant it is easy to 
see that we need H only up to the first order in e. Then there is only one term 
of V that is relevant: 

U(3) = (42) 

To compute the renormalization of the quadratic terms in the Hamiltonian we 
have to assign, according to the general procedure of Section 2, two A’s to be 
“high” and one “low”. Therefore, the relevant part of looks like Q 

^(3) ^ + d,Al^A\,Al^) (43) 

We now write A^ and E^ in second-quantized form, 

= E (44) 

= E*(t)" 

k 


where the creation and annihilation operators satisfy 




6t 

i "> ^pj 



kikj\ 
fc 2 ) 




(46) 


After normal-ordering, VA) can have three kind of terms: aj.a|,, OfeUp and aj^ap. 
As was explained in Section 2, only the hrst two types will be used to in order 

^We use the following convention: fields with subscript 1 indicate “low momenta” fields 
and fields with subscript 2 indicate “high momenta” fields 
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to determine so the final form of is therefore 
^(3) ^ ^fatc ^ f ^3^ 

P.9 

N ' j Zi-^ujp 

In order so solve equation i) with the interaction given by (^) we use 
H2 = Y^ uj{k)al\ai,, [i?2, OpI] = ujpa^l, [i?2, Qp,] = -Wpa“,. 


( 47 ) 


(48) 


Then, 


P,<1 


^QJ + ap,ag,e J A^(x)j ^ ■ 

According to the equation (|ll|), the next Vt will be 

=_j epbc ^ f d3cr{2(a;Ja^tg-z(p+,).+„a^^c^.g*(p+,).^ [H,,djAtix)] 

tq, [i7i, A^^(x)]} 


(49) 


(^iOp -\- ^g)^ 


(50) 


To study the renormalization of the metric we have to determine the correc¬ 
tions to the term in the Hamiltonian. There are only two possibilities. The 
first one is — ^ [[Hi, Hq] , Hi], since Hq ~ A and Hi ~ E. The other comes from 
the normal-ordering of the quadratic part of the a term in the equation (^^. 
When computing the double-commutator, there is only one divergent term, 

-i[[i7i,Ho],Hi] = e2^i[i7i,A,'(x)] [i7i,AS(y)]r'7'""'Ma X 

k,p 

{SkT^ps + &;^) PAk)PM (e-^('^+p)^-v) + E<^^+p^^-y)) x 

(51) 


4 (Wfc -I- UJp){uJq + UJs)'^^U;kU!pUlqUls 

At this point we can say that momenta {k -|-p) and —{k+p) are essentially the 
momenta of the “low” fields. Therefore \k + p\ < p. We can now change the 
summation by the following trick: say k + p = r and |r| < p. Then k = r — p 
and the summation goes over k and r. Any possible divergence can only come 
from the summation over k. Using J“'^'=ybdc _ ^jig expression (27) can be 

simplified as follows, 


E 

r,k 


CA [i7i,A“(r)][iJi,A“(-r)] 




{uJk + UJpY 


[kikoi kiPo,^ Pj-y(^k^Ej^(^p^ (^ 2 ) 
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Noticing that the leading divergence of the expression is logarithmic, which 
means that we can neglect the difference between k and p (for the divergent 
contribution only) and using 


Pj~^{k)Pj^{p) = "I j^2p2 ^ 2, 


(53) 

(54) 


we obtain 


CA 

24 


J2[H„Anr)] [H„AU-r)]Y,^ 

r h 


CA 1 , 

-24 2 ^^'°’ 



d^xE^ (55) 


As we mentioned earlier, there is another term that can correct the E^ term 
of the effective Hamiltonian. It is the “Coulomb interaction” term from the 
kinetic term in (|3^). Looking at the order correction part of the (^ and 
choosing the A’s to be “high” and the E’s to be “low” we have 


f / A^A^)Ai^(y)G(^,y)m^)J^^Ay) (^6) 

As written, this term is not normal-ordered with respect to the high-energy 
vacuum. Using (|4^),(^ and (^ ) we obtain 

P'^G-y)G{x, y)EUx)El^{y) (57) 

Upon using the definition (|^ of G{x, y) one can see that the leading divergence 
is logarithmic and that the final expression reads 




d^x 


(58) 


The sum of the terms (M) and (|58| ) gives the total correction to the kinetic 
energy at one loop: 


5He^ = J d^xEf{x) (59) 

At this point we can say that the factor is therefore, 

Since the operator Ei is represented by a variational derivative with respect to 
gauge field, one would naturally expect that Ze^ should be equal to the 
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Z ^2 • It is therefore an important check on our method to show that it is indeed 
so. To compute ^32 = Z 3 we need to find out correction to the effective 
H. In comparison with the computation of the E'^ correction it is much more 
involved due to the fact that for most of the terms comes as a sub-leading 
divergence. We will not present detailed computation but sketch the main steps 
and give the final result. B^ contributions can arise from the following terms: 
I [V, flo] and normal ordering of the a term from the equation again. This 
term is similar to the the only difference being how the “high” and “low” 
components are assigned: 

( 61 ) 

Leading divergence for both terms is quadratic and gives correction of the form 
. The appearance of this term is related to our choice of the cut-off procedure 
as a way of regulating the theory; it can be dealt with by introducing a 
counter-terms in the bare Hamiltonian and defining appropriate boundary con¬ 
ditions at the ends of the RG flow trajectory |l^. To capture the logarithmic 
contribution one has to expand the denominators of the Hq up to the second 
order in the momenta of the “low” fields. Tedious but straightforward compu¬ 
tation gives 


= (Ef)( 62 ) 

The logarithmic divergent part of the normal ordering of the two from ( |^ 
gives 

( 66 ) 

The final result is 

which makes the corresponding Z factor 

= (65) 

This coincides with the value of Z 3 for QCD in the Coulomb gauge obtained in 
the Lagrangian formalism . 


5 Three point function renormalization 

Let us show briefly the renormalization of the three point function within our 
formalism. To do that let us recall the general method of Section 2. We will 
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have to be a little more careful in the analysis of the relevant contributions to 
the renormalized Hamiltonian. 

According to the notation of Section 2, we write the Yang-Mills Hamiltonian 
in terms of the low momentum and high momentum fields as: 


H = Hi+H2 + {vP + vP + + • • •) ( 66 ) 

where Hi is the part of the Hamiltonian that only contains low momentum 
fields, H 2 is the part that only contains high momentum fields, and the Y’s 
are the “mixing” terms. For convenience we have separated these last terms 
according to the vertex number and the high momentum creation-annihilation 
operators structure: 


= 






vP = e 


UJq 


(67) 


P,<? 


(‘ 






( 68 ) 


2-y^CcJp UJq 
iip+q)x'^ 

+ + ayp-dP+'l')-^ (^fabefcde^jadejcbe^ {x) A'^ (x)} X 




pi qi 


P,g 


LOq 


(69) 


Now we have to analyze which terms contribute to the three point vertices. At 
this order only fl at order e is needed, moreover, one can convince oneself that 
only rig, i.e, the first iteration of the unitary transformation leads to divergent 
contributions. Then eq.(||) reads in this case: 

H' = Hi+H2 + VpPvP+--- + i[Hunii]+i[H2,no\ + 
t[vP,no] + t[vP,no] + tp^\no]-^[[Hi,no], Og] - 
\[[H2M ^ 0 ] - ^ i [[VP , Og], f]o] - 

i , f^g], Ho] - ^ [[[i7i, Ho], Ho], ^ [[[ 772 , Hg], Ho], Ho] + • • • (70) 

As explained in Section 2, we choose Hg in in such a way that its commutator 
with H 2 cancels the mixing terms that contain high momentum annihilation 
operators or high momentum creation operators, but not both. Also, we note 
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( 71 ) 


that up to this order, only is of order e. Then Qq satisfies: 

[H2,no] = ivf 

(which is precisely the equation that gives eq.(^) ). 

Now it is not difficult to individualize the only terms that contribute to the 
three point vertex: 

= AH,+ + i[vf,nof^ ( 72 ) 

Here AH, stands for the normal ordering contribution (tadpole diagram) of 
order in the kinetic energy. In fact, at this order, the kinetic energy has a 
term of the form: 


^3 £abc xade remn 

e3 = e / / / 


A>l{x)G{x,y)A^^{v)d,G{y,z)A^{z)Et{x)El{z) 


(73) 


which, when properly contracted, generates the contribution: 


1 1 

AH, = -CA e'^ log 
6 


j^d,A<;A\A‘^^. (74) 


The remaining terms in eq.(72), are computed similarly to the ones of Section 
2. After a very lengthy, but straightforward computation, we get the following 
results: 



Finally, adding up all the contributions we have 


e’ (^) /“■“ I a.A‘A‘Al ( 70 ) 

This result corresponds to a renormalization constant Z, equal to: 


^1 = 1- -^CA T^log ( — 


12 


27T^ 


m 


(77) 


which is the same as the Coulomb gauge result in the Lagrangian approach as 
well [ p^ . 
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6 QED 

In this section we will outline a similar computation for Quantum Electrody¬ 
namics. The QED Hamiltonian can be written as follows 

i/ = i y d^x -|- J d^x I'i/j • d + mj ip + e ip A ■ 7 i/;|. (78) 

The imposition of the Gauss law constraint, 

diE^ = eip'j^ip (79) 


generates an “instantaneous” Coulomb interaction and E[ takes the form 


H = lJd^a:{E^ + B^)+^ 


ip{x)j°ip{x)G{x, y)ip{y)'y°ip{y) 


X Jy 


-\- J d^x ^ip ■ d + rnj ip + eipA ■ ^ip^ 


(80) 


where G{x, y) is given by the (^). According to the general idea we are supposed 
to split it in to “high” and “low” energy parts, 

i7i = - y d^x -I- Bi) + y d^x ipi (^id ■ 7 -I- ipi + eipiA ■ "/ipi{81) 

H2 = 

k p 

= ej d^x { 1 P 2 A 1 ■ ^ip 2 + eip 2 A 2 ■ jipi + eipiA 2 ■ ^ip 2 } (82) 


= e'^i’2ix)"f°ipiix)G{x,y)ipi{y)j°ip2iy) 

J X J y 


(83) 


Here dots mean that we show only those terms that play role in one loop effects. 
Using arguments similar to those in Yang-Mills theory one can see that in order 
to determine Z factors up to the order we need H only up to the first order. It 
turns out that only first two iterations in y/A are needed - Hq and Hi. In parallel 
to the previous sections, we determine Z 3 and Zi by identifying corrections of 
the type E^, B^ and ipA-^ip. There is only one commutator that can contribute 
to the term: | [V, Hq] where one has to take subleading divergence to identify 
the log corrections. E^ correction is given by two commutators Hq], Hi] 

and commutator of the “Coulomb” term with Hq : The final result 

is 


SHb2 = J d^xBf (85) 
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which leads us to the well-known answer for the Za-factor in QED: 



( 86 ) 


Up to this point all our Z's were identical to those known from the covariant 
calculations in the Coulomb gauge. However, when computing the other two 
renormalization constants: Z\ and Zi we find results which are different from the 
covariant ones. In the case of Zi - fermion kinetic term renormalization there 
are three possible contributions arising from |[U, flo], — 5 and 

normal-ordering of the four-fermion term in the Hamiltonian (|80|). Extracting 
■07 • 9'0-type corrections from each of these terms one can see that they cancel. 
The similar thing happens in case of Z\ as well: all the 07 • H0-type terms 
cancel at the order. This makes both Z\ and Zi equal to 1 at one loop. 
Nevertheless there is no contradiction between our result and the conventional 
one. In covariant formalism there exists a Ward Identity Z\ = Zi which is 
essential for the maintaining the gauge invariance of the effective action. It is 
satisfied, presumably, in our case as well. But values of Z\ and Zi are gauge 
dependent and cancel out from the final expression that defines the beta function 
for QED, which is determined by Z-j, only. There is an obvious reason why, say 
Zi must be 1 in our case. Since 0 and 0 are conjugated variables, in the 
Hamiltonian formalism one should be represented by the variational derivative 
with respect to the other. Then, similarly to the and BZ- terms for the gauge 
field, they will have inverse renormalization factors which will cancel each other 
in the final expression for the kinetic term for fermions. 

7 Geometric interpretation of the asymptotic 
freedom. 

Let us recall the Yang-Mills Hamiltonian, after integrating the high momenta 
modes down to the scale /i (at order e^), incorporating the correct renormaliza¬ 
tion factors Zi. 






(87) 


Similar expression for the QED Hamiltonian will be: 



Here is a fixed quantity at certain scale and all the dependence on the scale 
is hidden in the Z factors. The definition of the “renormalized” fields through 
the incorporation of the Z factors was done in analogy with the Lagrangian 
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covariant approach where the renormalized quantities are included in such a 
way that the “renormalized” effective action gives finite results when the cut-off 
is removed. But in the Schrodinger picture this requirement is not necessary 
as the fields are only coordinates of the configuration space and do not enter 
explicitly in the computation of correlation functions. With this fact in mind 
we will alter this requisite and adapt it to our needs. 

The usual covariant renormalization program puts the emphasis on the inter¬ 
actions, and the scaling properties of the theory are extracted from the study of 
the /3-functions. However, in a Hamiltonian description the kinetic term plays a 
significant role since essentially it is nothing but the Laplacian in configuration 
space. Hence, many of the properties of the QFT can be inferred from the geo¬ 
metrical features of the configuration space (as compactness, boundness, etc.). 
In particular we are interested in the change of the configuration space metric 
under the renormalization group flow. We claim that in asymptotically free the¬ 
ories the distance between configurations increases as we move to the UV limit, 
thus “flattening” the potential energy and consequently fading the interaction. 
We will support our claim with the analysis of the one-loop Yang-Mills theory. 

It is clear then, that in the spirit of our work we want to stress the kinetic 
term (better, the configuration space metric) over the potential energy and try to 
understand the asymptotic behavior of the theory through the renormalization 
flow properties of the distance in the configuration space. For this reason we will 
rescale the fields in such a way to transfer the renormalization scaling properties 
to the kinetic term. 

Then let us rescale the fields as: 

^ , Et ^ efl Z-Hp)Z,{p) (89) 

so the QCD Hamiltonian takes the form: 



where we have used the Slavnov-Taylor identities (adapted for the Hamiltonian 
formalism) n 


Zi ^ ^4 

Z^ Z\ 


(91) 


Note that with this normalization we have “homogenized” the potential term 
(up to an overall factor) by transferring all the cut-off dependence to the kinetic 
term. Now we can read from the kinetic term, the cut-off dependence of the 

^ By writing both the Hamiltonian H and Gauss law Q°' in terms of renormalized fields 
with the explicit Z factors and requiring [H, = 0 one can get linear relations between Z’s 

that lead to (pl|). 
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(inverse of the) metric. In fact, using the result of Section 3, we can write the 
cut-off dependent configuration metric as: 


G{ia)Ub){x,y:fi) = + ^ + (92) 

where is the metric defined in equation (^^ written in terms of the 

rescaled fields ( ^^ and cr. 

Looking at the QED Hamiltonian and using Zi = 1 and Z 2 = 1 one can 
see that no rescaling is needed at all, since all the cut-off dependence is already 
shifted to the kinetic term. 


Hqed = ^ ((^) +ipii-dtp + criIjA --fil)--- (93) 


Following our procedure the corresponding metric will be 

G(ij){x,y,y) = ^1 - i ^ log(/r/mfl)^ S(^ij)S{x - y) + 0{e^R). 


(94) 


At this point it is easy to compare relative behavior of the two metrics under 
the renormalization flow. Equation (92) clearly shows that the distance between 
configurations decreases as the cut-off is lowered, while the corresponding ex¬ 
pression ( p^ increases, sustaining then our claim. Incidentally, it is worthwhile 

to mention that the combination ^ that appears in equation ( 1 ^) is precisely 


the one that defines the /3-function of the Yang-Mills theory and, at least in the 
Lagrangian approach, it has been proved to be independent of the gauge fixing 
condition. Similarly, QED expression involves only Z 3 , since Z\ = Z 2 is I to 
this order, and Z 3 is the only constant that determines beta-function for QED 
and it is known to be gauge-independent to all loops. 


8 Summary and Conclusion 

There are several issues that are more natural to address in the Hamiltonian 
picture than in the usual covariant Lagrangian formalism. One of them, which 
we are interested in, is the relevance of the geometry of the configuration space 
to the properties of the corresponding Quantum Field Theory. The reason is 
simple: in the Hamiltonian formalism the kinetic energy term is nothing but a 
Laplacian operator in the configuration space and its topological and geometrical 
features determine the nature of its spectrum. Then it is natural to ask what 
is the behavior of the configuration space as we integrate out high momentum 
degrees of freedom. 

To answer a small part of this question was the aim of this paper. To be 
precise we were interested in the following aspect of the problem: the evolution 
of the distance between field configurations (and more precisely the metric) 
with the renormalization group in asymptotically free theories. In particular 
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we state the following conjecture: in asymptotically free theories the effective 
distance between configurations decreases as high momenta degrees of freedom 
are integrated out. 

To support this statement we first developed an original renormalization 
group technique for Hamiltonian formalism in the framework of perturbation 
theory. This method resembles the Hamiltonian renormalization approaches of 
Glazek and Wilson and Wegner Q and operates by a progressive diagonal- 
ization of the Hamiltonian by means of a succession of iterative unitary transfor¬ 
mations followed by a projection onto the Hilbert space of the low-momentum 
degrees of freedom. Finally we applied the formalism to two conspicuous QFT’s: 
Quantum Electrodynamics and Yang-Mills theory in 3-1-1 dimensions, where we 
constructed the renormalized Hamiltonian up to one loop. 

Our results were substantially supportive of our conjecture. In the case of 
Yang-Mills, an asymptotically free theory, the one-loop metric renormalization 
showed that in fact the distance between configuration increases as the momen¬ 
tum scale increases, and on the contrary for QED, not asymptotically free, the 
behavior of the metric is the opposite. 

We are aware, of course, that our results are not decisive but just consistent 
with the conjecture. After all we have only studied two examples at one-loop 
order in perturbation theory. However from the examples considered we can 
observe a pattern that seems to repeat at any instance: when moving all the 
weight of the renormalization group onto the configuration space metric, it 
acquires a renormalization factor which is a function of the same combination 
of renormalization constants that defines the /9-function of the theory, and thus, 
presumably, inheriting its asymptotic behavior properties. 


Acknowledgements 

We want to thank Prof. V.P. Nair for suggesting this problem and his advice, 
encouragement and input during the work on this project. We are also grateful 
to Prof. P. Orland for helpful discussions and critical reading of the manuscript. 
G.A. was partially supported by CUNY RF grant 6684591433. E.F.M. was 
supported by CONICET and CUNY Collaborative Incentive Grant 991999. 


Appendix A 

In this appendix we show that the ambiguity in the definition of the effective 
Hamiltonian due to the freedom in the solution of equation (H) (and generally 
is just the standard ambiguity of the Hamiltonian operator, namely the 
freedom of unitary transformations. 

Let us recall that the solution of & or (|^) is not uniquely defined; if 
Qn is a solution so will be + On with On satisfying the homogeneous equation 
[H 2 ,On] = 0. Now let us consider two sets of H’s - say and fln\ all of 
them satisfying the proper equations (||) and ( |l^ but generated from different 
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type of solutions. Then we have 


Ha 

Hb 


^ — iO, 


(a) 


—iO., 


-iO.', 


... e 


He 


in: 


(a) 




(a) 


ini'‘'> An 


(t) 


(95) 


Obviously Ha and Hb are unitarily related: Ha = U~^HbU and con¬ 
sequently they have identical spectrum. What we have to show is that this 
property remains a/ter projecting onto the “high” perturbative vacuum \0high)- 
Following the prescription given in section 2, we have shown that, up to 
given order n in coupling constant A and a given order m in A we can write: 

Ha = i?l“^(low) + i7fr-(high)+ ^ alsi;^ap + 0{X-,{^l/AD 

k,p>^ 

Hb = i?f>(low) + i/fr-(high)+ ^ 4^g^ap + 0(A^(M/A)™) (96) 

k,p>g 


where the Hamiltonians H^'^^ depend only on the low energy modes, H 2 is the 
free Hamiltonian for the energy modes and are arbitrary operators of low 

and high frequency modes. 

Without losing any generality we can assume that h[°'^ ) and H^’’^ are diag¬ 
onal, as they can always be brought to that form with a “low energy modes”- 
unitary transformation that respects the structure of (p^. 

Now we will show that the low energy Hamiltonians h[°'^ and have the 
same spectrum and consequently they are unitarily related. 

Consider the eigenvalues equations for h[°'^ and h[^'^: 

c = K C , h[^^ = Ei (97) 


Using standard perturbation theory we can compute the eigenvalues of the 
whole Hamiltonians, A and B as an expansion in powers of the matrix elements 
of the interaction terms 

^ (98) 

k,p>fi 


We get, for the eigenvalues of the operator Ha, 

'a, n|U“| 7 , to) ( 7 , to|U“|q;, n) 


=EDEDMVDn 


-E- 


+ E° - E°‘ - E° 

■^'v ' -^m. -^71. 


(99) 


and a similar equation is valid for the eigenvalues of the operator Hb ■ But the 
low energy spectrum corresponds to those states with n = 0, and in this case, 
due to the particular form of the interaction, all the perturbative contributions 
vanish and the eigenvalue E^ is the exact eigenvalue of the whole Hamiltonian 
(up to the given order in A and /i/A): 

El^o^ = E^^ + o D, DAD ( 100 ) 
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And finally, since both Hamiltonians Ha and Hb have the same spectrum, 
or since is the same for both Ha and Hb, we deduce that 


E 


a 

<y 


A^ + 0(A-,WAn . 


( 101 ) 
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